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We study photon-assisted tunneling (PAT) through weakly-coupled doped 
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J^ , semiconductor superlattices (SL's) in the high voltage regime. A self- 
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f^ , consistent model which treats the Coulomb interaction in a mean field ap- 

proximation is considered. We discuss the formation of electric field domains 
in the presence of THz radiation and the appearance of new multist ability 
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multiple photon emission and absorption. 
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Recently, the transport properties of semiconductor nanostructures in the presence of AC 
fields have been the subject of interest ||l]-§[. Among the experiments based on the response 
of a nanostructure to an AC field one might consider PAT, electron pumps, turnstiles, 
and others . On the other hand the transport properties having their origin in Coulomb 
interaction have attracted a great deal of attention in recent times. In weakly-coupled SL's, 
multistability, electric field domain formation, self-sustained current oscillations and chaos 
have been much studied |@,||. In this work we deal with PAT in weakly-coupled doped SL's 
whose transport mechanism is sequential tunneling. In this regime the intrinsic miniband 
width is typically one order of magnitude smaller than the scattering induced broadening 
(for the samples we are going to analyze the miniband width is ~ O.lmeV). In this case the 
miniband transport can be neglected and the transport at low bias voltages is governed by 
sequential resonant tunneling from ground state to ground state in adjacent wells. The PAT 
through a SL in the low bias regime, where the non linear effects of the Coulomb interaction 
are small, has been explored both experimental |I|] and theoretically ||^. In this regime, 
and for certain values of the intensity and frequency of the AC field, the sample can exhibit 
negative conductance []I],^. Here, we concentrate in the high bias regime where the Coulomb 
interaction has to be considered in order to discuss electric field domains formation in the 
presence of THz fields. We study the formation of new multistability regions induced by 
the AC field and we show how the electric field domains can be supported by the virtual 
PAT sidebands 0. An extension of the model of Ref. 10 is put forward in order to account 
for the effect of an AC field on the I-V curve of a doped SL. The model assumes that the 
characteristic time of intersubband relaxation is shorter than the tunneling time and then 
that only the ground state of each quantum well is populated. In order to model relaxation 
in the transport direction due to scattering , we assume that the spectral functions in the 
wells are Lorentzians. Relaxation in the planes perpendicular to the transport is taken into 
account by imposing current conservation through the whole structure that will determine 
the sequential current as well as the Fermi energies within each well. The assumption of 
having Fermi distribution functions in the wells implies some scattering mechanism that 



thermalizes the electrons towards an equilibrium state P]. If there is no mechanism for 

obtaining thermal equilibrium the nonequilibrium transport must be treated with a suitable 

quantum kinetic technique (for example the Keldysh technique [Q). 

The currents are calculated by means of the Transfer Hamiltonian method (TH) that consists 

in connecting by time dependent perturbation theory the different parts of the structure 

(decoupled in the remote past) by tunneling matrix elements p. 

The effect of the AC field is included in the single particle energies in the different isolated 

regions of the structure (leads and wells); ei(t) = ej + eVi(t) before the tunneling couplings 

are switched on [Q. Due to this time dependence the single particle propagators acquire 

phase factors. The noninteracting retarded Green's function at the i-th region becomes: 
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In this expression J„ is the Bessel function of first kind and the time dependent voltage is 
Vi{t) = Vi cos tut. Following the TH method we obtain the transmission probability from the 
z-th well to the i + 1-th well , from the emitter to the first well and from the A^-th well (A^ 
is the number of wells) to the collector. From the transmission probability we evaluate the 
time averaged sequential current . The interwell current is Q: 
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where e^^j- is the j-th resonant state of the i-th well {rimax is the number of subbands partici- 
pating in the transport) , Tj(e, e + mhuj) is the inelastic transmission through the i-th barrier. 
In the argument of the Bessel functions , Fd = Vi^i — V^ is the potential drop between the 
i-th well and the i + 1-th well due to the time dependent field; F the intensity of the time 
dependent external field and d the period of the SL. The current from the emitter to the 
first well, Jo,i, and the one from the A^-th well to the collector, Jn,n+i, are also derived in 



our model. The electrons in each well are in local equilibrium with Fermi energies e^j^ which 
define the charge densities rij. For a given set of variables {e^^} and in the stationary regime 



the currents have to fulfill the set of equations |Tn| 



Ji_i,i - Ji,i+i = i = l,...,N. (3) 

Here Jj,j+i = Ji,i+i(e^^, e^^^,, $), Jo,i = JoA^uji,^), and Jn,n+i = Jn,n+i{(^ujn,^)- $ denotes 
the set of non linear voltage drops (potential drops at the accumulation and depletion layers, 
barriers and wells) through the structure caused by the accumulation of charge densities rii 
in the wells. These potential drops as well as all the quantities in the problem are calcu- 
lated self-consistently for each applied voltage including the Coulomb interaction in a mean 
field approximation ||10[. Our numerical procedure allows us to obtain both the stable and 
unstable solutions as well as all the multistability regions in the I-V curve. 
In Fig. 1 we plot the I-V curve of a SL consisting in 10 wells with 90AGaAs wells and 
40AGao.5Alo.5As barriers. The doping at the leads is N^ = 2 x lO^^cm"'^ and in the wells it 
is N^ = 1.5 X lO^^cm"^, the half-width of the resonant states is 7 = 2meV and T=0. 
After the initial peak which is determined by Cl-Cl sequential tunneling {Cj are the reso- 
nant states of the wells ordered starting from the ground state) through the whole SL the 
current evolves along a series of branches (solid lines). This behaviour is explained by the 
formation of a charge accumulation layer in one of the wells (domain wall) [^,^ that splits 
the SL in two regions with low and high electric field respectively (see inset in Fig. 2). In- 
creasing the voltage, this charge cannot move continuously through the SL. This motion can 
only occur for voltages allowing resonant interwell tunneling, in this situation the domain 
wall moves from the i-th well to the i — 1-th well. In the I-V characteristics this process 
leads to a series of sharp discontinuities (the stable branches are connected by unstable ones, 
dotted lines in the figure). In the static case, transport in the high electric field domain is 
only possible by C1-C2 resonances, whereas in the case with AC the resonances occur be- 
tween the virtual PAT emission sidebands corresponding to the C2 subband in the i + 1-th 
well and the PAT absorption sidebands corresponding to the CI subband in the i-th well . 



These processes are shown in Fig. 2 where the I-V curve in the presence of THz radiation 
{F = 0.47 X lO^V/m, uj = 3THz) is plotted. In the inset we have depicted the calculated 
potential profile of the SL for a voltage Vq = 0.86V. The inset shows that the virtual PAT 
channels can support electric field domains (the dashed lines represent the PAT channels), 
the blow up shows that the high field domain is supported by Cf-C2 tunneling involving 
absorption of two photons. Increasing the intensity of the AC field to F = 0.95 x lO^V/m 
(Fig. 3) the probability of having multiphotonic effects increase, leading to multistability 
of the branches. The inset shows a magnification of the first branch, the circles mark the 
stable operating points for a fixed voltage. At Vi = 0.16V transport in the high field domain 
occurs via tunneling between the two-photon absorption virtual state associated with CI 
and the two-photon emission virtual state associated with C2. At V2 = 0.19V the branch 
develops a multistable solution (five solution coexist, three stable, two unstable). These 
solutions correspond to a different number of photons emitted in C2: one photon in the 
highest current stable solution (circle a), two photons in the lowest current stable solution 
(circle c); the process from the highest current to the lowest one involves the motion of 
the domain wall. The situation repeats periodically as the domain wall moves, giving the 
sawtooth structure in the current. 

In summary, we have presented and solved a microscopic self-consistent model for the se- 
quential current through a weakly-coupled doped SL in the presence of THz radiation. We 
have shown how the electric field domains can be supported by the virtual PAT sidebands. 
For strong THz fields multiphotonic effects lead to appearance of new multistability regions 
in the I-V curve. 
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FIGURES 
FIG. 1. I-V characteristic for a SL without AC. The continuous (dotted) hues correspond to 

stable (unstable) solution branches. 

FIG. 2. I-V characteristic for a SL with AC (F = 0.42 x lOf'V/m, uj = 3TiJz).The inset shows 
the calculated potential profile at V=0.86 V. 

FIG. 3. I-V characteristic for a SL with AC {F = 0.95 x IQ^V/m, uj = 3riIz).The inset shows 
a blow up of the first branch. 
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